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A statistical model to explain the gamma-ray variability and flares of the Crab nebula 
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Abstract: Recently the AGILE and Fermi/LAT detectors uncovered giant 7-ray flares from the Crab nebula. The duration 
of these flares is a few days. The Fermi/LAT data with monthly time binning further showed significant variability of 
the synchrotron tail of the emission, while the inverse Compton component was stable. The simultaneous or follow-up 
observations in X-ray, optical, infrared and radio bands did not find significant flux variation. Based on these observations, 
we propose that the 7-ray variability and flares are due to statistical fluctuations of knots that can accelerate electrons to 
~PeV energies. The maximum achievable energy of electrons is adopted to be proportional to the size of the knot, which 
is assumed to follow a power-law distribution. Thus the low energy electron flux will be stable due to the large number 
of small knots, while the high energy electron flux may experience large fluctuations. Monte Carlo realization of such a 
picture can reproduce the observational data quite well given proper model parameters. 
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1 Introduction 

The Crab nebula is one of the most famous source in 
the sky, and is the most typical multi-band laboratory. It 
has been widely studied in all wavelengths, from radio 
to very high energy (VHE) 7-rays. The broadband non- 
thermal emission can be well modeled with a leptonic sce- 
nario, with the synchrotron radiation from the radio to GeV 
7-rays and inverse Compton (IC) radiation above ^GeV 
Q] . The transition from synchrotron to IC component was 
clearly seen by CGRO/EGRET @ and Fermi/LAT 0. 
Two populations of electrons, i.e., the radio electrons and 
wind electrons are required to fit the data HI . The overall 
emission from the nebula seems to be stable in radio, op- 
tical, X-ray and VHE 7-ray bands, thus it was thought to 
be a "standard candle" and was often used to calibrate the 
detectors. 

However, years ago people began to find that the MeV- 
GeV emission of Crab nebula was actually variable instead 
of stable, from the COMPTEL and EGRET observations 
ID 13. Furthermore the detailed images of the Crab neb- 
ula in optical and X-ray bands also indicated dynamical 
structures at small scales. The HST observation revealed 
highly variable wisps and knots in the nebula |6|. X-ray 
observations by ROSAT and Chandra uncovered a jet-torus 



structure of the inner nebula, which was also dynamically 
variable (5]Q. Recently in September 2010, the AGILE 
X-ray satellite detected a 7-ray flare above 100 MeV from 
Crab nebula, with a flux enhancement of ~ 3 times higher 
than the average one and a duration of about 3 days [8]. 
This flare was confirmed by the Fermi/LAT detector, with 
an even higher flux enhancement [9]. The search for the 
archive data of AGILE and Fermi/LAT further unveiled 
other flares, in October 2007 for AGILE and in Febru- 
ary 2009 for Fermi/LAT [8J |9), which may indicate that 
such flare events occur with a timescale about one year. 
Most recently in April 201 1, AGILE detector observed an- 
other super-flare from Crab nebula iflOl . Furthermore the 
monthly binned light curve of Fermi/LAT data showed that 
the synchrotron component is variable, however, no varia- 
tion was found for the high energy IC component |9 1. 

There were many other simultaneous or follow-up observa- 
tions for the flare in September 2010 in other wavelength 
bands. In the VHE 7-ray energies, ARGO-YBJ collabo- 
ration reported a detection of a flux enhancement around 
TeV with a possibly longer duration (ATel #2921). The re- 
sults from MAGIC and VERITAS telescopes with limited 
exposure time, however, did not find any flux variation dur- 
ing the flare phase (ATel #2967, #2968). In X-ray, optical, 
infrared and radio bands, no significant flux enhancement 
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was discovered either (ATel #2856, #2858, #2866, #2867, 
#2868, #2872, #2882, #2889, #2893, #2903). 

Some theoretical models are proposed to explain the flare 
event lfTTlfT2l . In lfl3l we proposed that the 7-ray variabil- 
ity and flare were due to the statistical fluctuation of the 
acceleration units which were responsible for the highest 
energy electrons. It is natural to expect that events that can 
generate the highest energy electrons are rarer, and there- 
fore would suffer from the largest fluctuations. The lower 
energy electron spectrum should be steady-like because 
many more accelerators can contribute to them simultane- 
ously. Since the synchrotron 7-rays are produced by the 
highest energy electrons (~PeV), while the synchrotron 
emission at lower energies and IC emission at higher ener- 
gies are produced by lower energy electrons, we can easily 
explain why we only detect variability in Me V-GeV 7-rays, 
neither in lower nor in higher bands. 

2 Model 

We assume that electrons are accelerated in a series of 
knots, which may relate with the magnetic hydro-dynamic 
turbulence of the plasma. For each knot the output energy 
spectrum of electrons is a power-law function Fi(E) oc 
E~ ac with a maximum electron energy -E* nax , which is 
proportional to the size of the ith knot. It is further as- 
sumed that the size of knots has a power-law distribution 
P(ri) oc r^* 3 . The total electron spectrum can be got 
through adding the contribution from all the knots together. 

The comoving system synchrotron spectrum from a knot 
with size is v'F v i oc v'~ a " cxp(— v' /v' max ), where 
v' max oc (i?n iax ) 2 oc rf, and a v = (a e — 3)/2 is the syn- 
chrotron spectral index, with a e being the power-law in- 
dex of the electron distribution. Changing to the observer's 
frame, there will be a frequency shift v' — > 5v' and a flux 
enhancement F v i — > 5 3 F V >, where 5 is the Doppler factor 
of the knot. 

It was shown that the maximal energy of synchrotron emis- 
sion in the magnetic field dominated acceleration regime 
is ~ 70 MeV, which is due to the fact that the electron 
may lose most of its energy in one cycle of the Larmor mo- 
tion [14]. For smaller accelerators, the maximum energy 
is further limited by the size of the accelerator. The ob- 
servational flare of Crab nebula has energies higher than 
GeV, which means that either there is Doppler boost of the 
flare event, or the acceleration of the electrons is not the 
shock-like scenario [9|. In order to explain the high en- 
ergy photons of the flare event, we therefore employ a mild 
Doppler boost factor. Only knots with large enough size 
and large enough Doppler factor can contribute to high en- 
ergy synchrotron radiation to explain the flares observed by 
Fermi/LAT. 

The cooling of the electrons need to be considered. The 
cooling time of the electrons is energy dependent. For syn- 
chrotron photon with energy e the cooling time is t c ~ 
1.5( J B/mG)- 1 - 5 (e/keV)-°- 5 J-°- 5 yr. If we adopt an av- 



erage magnetic field of the Crab nebula B w 0.1 mG, a 
mildly Doppler factor i5 ~1, the electrons correspond to 
e ~ 2.5 eV photons will have a cooling time comparable 
to the age of the nubela i agc w 10 3 yr. If the production 
of the knots is continous, the equilibrium electron spectral 
index should beQ a e — a l p + 1, for electrons with cool- 
ing time shorter than i ago . For the electrons with lower 
energies, whose cooling time is even longer than the age of 
the nebula, there will be no cooling at all and the injection 
spectrum will keep unchanged. Here we will focus on the 
high energy part, e.g., the synchrotron spectra from X-ray 
to 7-ray band, we will take a cooled spectrum of electrons 
as input of the model. 

3 Monte Carlo simulation 

We realize the above picture by a Monte-Carlo simula- 
tion. The injection spectrum of electrons from the knots 
is adopted as ajfj = 1.6, which corresponds to the fit to 
the radio-optical spectrum of the nebula fl5l . After tak- 
ing into account the cooling effect the electron spectrum 
is a e — 2.6. We directly generate the synchrotron spec- 
trum instead of the electron spectrum. The corresponding 
synchrotron spectral index is then a v = —0.2. The maxi- 
mum energy of the synchrotron photon in the knot comov- 
ing system is proportional to r\. We normalize i/ max of 
the largest knot(s) to be 55 MeV to account for the con- 
straints of cooling effect during the acceleration. The nor- 
malization of the synchrotron emmisivity from each knot 
is adopted a volume-proportional factor r 3 . As for the 
Lorentz factors of the knots, we assume a Gaussian dis- 
tribution with the mean value V = 2.0 and a standard de- 
viation a = 0.25. Such a Lorentz factor is consistent with 
the upper limit of the typical velocity of the jet |9|. The 
angle 9 between the knot motion and the line-of-sight is as- 
sumed to be randomly distributed. For 3<r range of T we 
find the Doppler factor S = 1/T(1 — v cos 9) is in the range 
[0.18, 5.5]. Finally the power-law index of knot size distri- 
bution is adopted as (3 w 4.8. After adding the contribution 
from all the knots with different sizes together we get the 
total synchrotron spectrum v' F v i oc v'~ 0/2 , which is can 
well reproduce the observed spectra of Crab nebula from 
optical to MeV 7-ray band, as shown in Figure Q] 

The simulated spectra together with the observational data 
are shown in Figure Q] For comparison we show the fit- 
ting results to the broad band data with the thick blue line, 
with a broken power-law electron spectrum. The elec- 
tron spectral indices are /3i//?2 = 1.60/3.45 below/above 
the break Lorentz factor 7b r = 1.1 x 10 6 . The high 
energy cutoff is adopted as a super-exponential behavior 
7-* exp[-( 7 / 7cut ) 5 ], with 7cut = 5 x 10 9 and 5 = 2.0. 
The magnetic field B = 124 /iG, which is a constant in 
the all nebula. This result is similar to the model invok- 
ing two population electrons as introduced in [ 1 1 and fl5l . 

1. Note that here the equilibrium spectrum actually represents 
a series of knots with the same sizes, instead of a single one. 
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Figure 1: Spectral energy distribution of the Crab neb- 
ula. Thin-green-solid lines: simulated synchrotron spec- 
tra in 30 realizations; thick-blue-solid: fit of the average 
multi-wavelength emission; black-dashed: examples of un- 
usual events which may be responsible for the flares. Data 
are from: radio-optical ED, INTEGRAL E3, COMPTEL 
E), Fermi/LAT 0, MAGIC EH, HEGRA ED, HESS 
1201 . Tibet- AS7 ED, and Fermi/LAT flares I9ll22l. 

The average of the simulated spectra is similar to the steady 
state fit, and is consistent with the long term Fermi/LAT 
and COMPTEL data. However, there is scattering of the 
simulated spectra, which reflects the observational variabil- 
ity of the 7-ray spectrum. In the X-ray band, the fluctuation 
is very small, which is close to a steady state emission. 

In this model, the IC component should not vary signifi- 
cantly. This is because the IC photons around TeV energies 
are primarily produced by the low energy electrons. For in- 
stance for ~ 100 TeV electrons the synchrotron photon en- 
ergy is typically 0.016(5/mG)(£/TeV) 2 keV - 16 keV 
for B ~ 0.1 mG, which suffers from very small fluctuation 
as shown in FigureQ] The corresponding IC photon energy 
is eic ~ max(7 2 e, 7m e ). For background photon energy 
e ^eV the IC photon energy is as high as 100 TeV. Even 
for the cosmic microwave background photon e ~ 10 -3 
eV we still have eic ~ 10 TeV. Therefore the fluctuation 
of the IC component should be indeed small. This result is 
consistent with the Fermi/LAT observations [9|. 

The two large flares in September 2010 and April 201 1 are 
not well reproduced in the simulation. We show two illus- 
trations of these events with proper parameters in Figure 
Q] For the September 2010 flare the model parameters are 
5 = 5.5, v' max — 70 MeV These parameters are reasonable 
in the present frame. However, for the April 201 1 flare, we 
may need 5 = 8.0, v' max = 70 MeV, which seem to be very 
extreme. Such a flare can be regarded as an event with very 
small probability. 

To investigate the fluctuations of emission in detail, we cal- 
culate the light curves of the simulated photon fluxes. Since 
there is purely statistical fluctuation of the fluxes, we can 
adopt independent realizations to represent each time bin. 
The length of the time bin is proportional to the number of 



Figure 2: Light curves of the simulated synchrotron flux 
above 100 MeV. The upper panel shows the results with 
weekly bin and the bottom panel for monthly bin. The av- 
erage flux of the realizations is normalized to the observa- 
tional one ~ 6.2 x 10 -7 cm -2 s _1 (9), as shown with the 
horizontal lines. Also plotted in the lower panel are the 
Fermi/LAT observational data with monthly bin. 

knots. We need a normalization of the absolute time scale, 
which is adopted that in one year there is an enhancement 
in flux by at least a factor of 5-6 for weekly bins, which 
may be responsible for a flare. We plot the weekly and 
monthly bin lightcurves in the upper and lower panels of 
Figure |2 respectively. The Fermi/LAT observed monthly 
light curve [9 | is also plotted in the lower panel of Figure|2] 
for comparison. It can be seen that the predicted monthly 
variability scale of the simulated results is very similar to 
that observed by Fermi/LAT. 

Figure [3] shows the histogram of the flux distribution. For 
the monthly bin results, we can see clearly the similarity 
between the simulation and the observational data. In the 
right panel of Figure [3] the histogram of flux distribution 
for weekly bin is shown. We can see that there is a signifi- 
cant concentration toward low-flux, which is distinct from 
a symmetric distribution around the average value. The de- 
tailed analysis of Fermi/LAT data can test such a predic- 
tion. 

4 Conclusion 

In this work we propose that the fluctuations of the elec- 
tron spectra at the highest energy end (^PeV) are responsi- 
ble for the variability of the synchrotron tail in GeV 7-ray 
band. The electrons are thought to be accelerated in a se- 
ries of knots, with a size distribution P(Vi) oc rT^ and a 
distribution of the Doppler factor. The maximal achievable 
energy of the electrons in the co-moving frame is assumed 
to be proportional to the size of the knots. Thus the very 
high energy electrons are generated by the very rare knots 
with both large sizes and high Doppler boosts, and hence, 
suffer from large fluctuations. On the other hand, the low 
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Figure 3: Histograms of the synchrotron flux distribution above 100 MeV, for monthly bin (left) and weekly bin (right), 
respectively. In the left panel the Fermi/LAT observation result is also shown (dashed histogram). 



energy electrons can be accelerated by many smaller knots, 
and the fluctuations are smoothed out. Therefore only the 
variability in MeV-GeV 7-ray band is observed, because it 
is the direct reflection of the highest energy electrons. The 
lower energy synchrotron component and the higher energy 
IC component are relatively stable due to the less fluctua- 
tion of the low energy electrons. The expected variability 
of the monthly bin fluxes above 100 MeV are well consis- 
tent with that observed by Fermi/LAT. The two large 7-ray 
flares can also be naturally accounted for without additional 
assumptions. 

It has been revealed recently that in the X-ray band the Crab 
nebula is actually no longer "standard candle", but experi- 
ences quasi-periodic variability with a level of several per- 
centage 11251 . Such large scale variability may be due to 
the overall energy injection of the nebula or the inhomoge- 
neous plasma flow induced magnetosonic waves [25]. The 
model proposed here may not be responsible for such vari- 
ability. However, it would be interesting to investigate the 
statistical fluctuations of the X-ray fluxes in different en- 
ergy ranges after removing the large scale variability, and 
compare with our model prediction. 

One thing we should keep in mind is that in the present 
model we do not expect the variability of the IC compo- 
nent. This conclusion is consistent with the Fermi/LAT 
data, however, in a relatively long time scale (monthly bin). 
During the flare phase, we note that the ARGO-YBJ ex- 
periment detected a potential flux enhancement in TeV en- 
ergy range (ATel #2921). Although it was not confirmed by 
the Cherenkov telescopes MAGIC (ATel #2967) and VER- 
ITAS (ATel #2968), we should note that the exposure time 
of MAGIC (1 hr) and VERITAS (2 hr) is very limited. It 
is thus very important to further monitor the Crab nebula 
in VHE energy band by e.g., Tibet-MD El or LHAASO 
E4l . to search for the variability, which can be a crucial 
test of the current model. 
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